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Abstract In vacuum high-intensity lasers can cause 
photon-photon interaction via the process of virtual vac- 
uum polarization which may be measured by the phase 
velocity shift of photons across intense fields. In the 
optical frequency domain, the photon-photon interac- 
tion is polarization-mediated described by the Euler- 
Heisenberg effective action. This theory predicts the vac- 
uum birefringence or polarization dependence of the phase 
velocity shift arising from nonlinear properties in quan- 
tum electrodynamics (QED). We suggest a method to 
measure the vacuum birefringence under intense optical 
laser fields based on the absolute phase velocity shift by 
phase-contrast Fourier imaging. The method may serve 
for observing effects even beyond the QED vacuum po- 
larization. 



1 Introduction 

To observe nonlinear responses of matter, the pump- 
probe technique is widely used: Matter is first excited 
by an intense laser pulse and then probed by a delayed 
weaker laser pulse. When the vacuum is considered as 
a part of matter, the most natural approach to probe 
it is, hence, the pump-probe technique. Maxwell's equa- 
tions in vacuum, however, allow only for linear superpo- 
sitions of laser fields. In quantum mechanics, a photon 
can be resolved into a pair of virtual fermions over a 
short time via the uncertainty principle in the higher fre- 
quency domain even below the fermion mass scale. The 
loop of the virtual pair provides a coupling to photons, 
resulting in a photon-photon interaction. In the optical 
frequency domain, the electron-positron loop and pos- 
sibly the lightest quark-antiquark loop are expected to 
give rise to the photon-photon interaction with the mass 
scale of the electron being 0.5 MeV/c^ and of the lightest 
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Fig. 1 Polarization directions of the target and probe lasers. 

quark ranging from ~ 1 — 100 MeV/c^, respectively. Be- 
low the electron mass scale, there is no known mass scale 
relevant for photon-photon interactions in the standard 
model of particle physics. In this paper we focus on the 
photon-photon interaction in the optical laser frequency 
range based on quantum electrodynamics (QED). 

In the low- frequency collision htu <C nieC^, it is suf- 
ficient to describe the photon-photon interaction by the 
effective one-loop Lagrangian [Tl[51|3] 

L,-ioo, = i^—MF,^Fn^ + nF,,Fn\ (1) 
360 m| 

2 

where a = j- is the fine structure constant, TOg is the 
electron mass, F^i, — dA^/dx^ — dA^/dx'^ is the an- 
tisymmetric field strength tensor and its dual tensor 
Pf"" = l/2e^""-PF,p with the Levi-Civita symbol e^""'P. 

Based on this Lagrangian, the dispersion relation for 
photons in vacuum is expected to be modified by intense 
electromagnetic fields. This effect under a constant elec- 
tromagnetic field was first discussed by Toll [T . At opti- 
cal frequencies, we may approximate the time-varying 
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electromagnetic field as a constant field, because the 
relevant time scale for the creation of virtual electron- 
positron pairs is much shorter than that of the inverse 
of optical frequencies. We can discuss the dispersion re- 
lation and the birefringent nature via measurements of 
the refractive index, i.e., the inverse of the phase velocity 
as illustrated in FiglU where a linearly polarized probe 
laser beam crosses a linearly polarized target laser beam. 
The measurements of the phase velocity shift when the 
electric fields of both lasers are either parallel or normal 
to each other are specified with subscriptions || or _L, re- 
spectively. The theoretical derivation of these quantities 
in the linearly polarized electromagnetic field of the tar- 
get (the so-called crossed-field configuration, where the 
electric field E and magnetic field B are normal at the 
same strength) was originally studied in [5J|S] and fur- 
ther derived from the generalized prescription based on 
the polarization tensor, applicable to arbitrary external 
fields, in [7^. This results in 

/ 1 8 2Pt 

.,/c = l-H,2^, (2) 
45 pc 

where W||/c and vj_/c are the phase velocities when the 
combination of linear polarizations of the probe and tar- 
get lasers is parallel and normal, respectively. The quan- 
tity Pc = ■mle'/fi^ - 1.42 x lO^J/pm^ is the Compton 
energy density of an electron and pt is defined as Zk / 
where k is the wave number of the probe electromagnetic 
field with the unit vector of k. The Lorentz-invariant 
quantity Zk is defined as 

Zk = (kc^F-nihF^), (3) 

and the relation to the energy density in the crossed 
field condition is 

^ = e\l + Ck-h)r, (4) 

with e = E = cB and h — B x E with " indicating the 
unit vector. Thus the second terms in Eq. ^ show that 
the deviation of the phase velocities of light iiy and v± 
are proportional to the field energy density normalized 
to the Compton energy density of an electron. The shift 
of the refractive index from that of the normal vacuum 
is on the order of 10^^^ for the energy density of 
1 3 / fim!^ corresponding to the power density of a high- 
power laser beam focused to lO^^W/cm^ at its waist. 
The refractive medium exhibits a polarization depen- 
dence, i.e., it shows birefringence. The difference in wy 
and v± in Eq. ^ results from the first and second terms 
in the bracket of the effective one-loop Lagrangian in 
Eq. ((H). 

The dispersion relation and the birefringence under a 
constant electromagnetic field in the UV limit (lu — > cxi) 
may be evaluated via the Kramers-Kronig dispersion re- 
lation, as discussed in [S] . The phase velocity in both UV 



and IR is expected to be subluminal (vphase < c) under 
the infiuence of the QED field [MZ]. The UV limit of 
the phase velocity is supposed to govern causality which 
should not exceed the velocity of light in vacuum. There- 
fore, it can be a fundamental test of a variety of effective 
field theories in the IR by testing whether the phase ve- 
locity in the UV limit, extrapolated from that of the 
IR, is superluminal {vphaseioo) > c) or not. Thus far the 
dispersion relation from IR to UV is theoretically known 
only in the QED field [S]. However, there is no data so 
far even in the domain of IR frequencies. It is impor- 
tant, therefore, for experiments to quantitatively verify 
or disprove the QED prediction. We note that the mea- 
surement of the refractive index in the domain of higher 
frequencies may be sensitive to the part of the anomalous 
dispersion where the real part of the refractive index rises 
as discussed in and, also, the measurement of the 
electron-positron pair creation [10.11. 12. 13 in strong 
electromagnetic fields may be directly sensitive to the 
absorptive or imaginary part. The Kramers-Kronig rela- 
tion connects the real and imaginary parts of the forward 
scattering amplitude or the refractive index. Therefore, 
the systematic measurements of real and imaginary parts 
over a wide frequency range may provide a test ground 
of QED and the Kramers-Kronig relation itself, when it 
is applied to the vacuum. 

The key issue is how to detect the extremely small re- 
fractive index change, resulting from the photon-photon 
interaction between the target and probe lasers. The con- 
ventional way in the X-ray frequency range is based on a 
measurement of the ellipticity caused by the target field- 
induced birefringence with respect to the linearly polar- 
ized incident probe photons [9ll23j . Since nowadays high- 
precision X-ray polarimetery technique is available |14| . 
we may reach the sensitivity to QED-induced birefrin- 
gence, if high-intensity lasers such as those attainable in 
ELI [TS] are provided. As explained above, the probe fre- 
quency dependence of the birefringence is important to 
complete the QED-induced dispersion relation. There- 
fore, we need measurements in the optical frequency 
range as well. The conventional ways in the range of 
optical frequencies that were performed [21j and pro- 
posed [22] are again based on a measurement of the el- 
lipsoid caused by the birefringence and a measurement of 
the rotation angle of a linearly polarized probe laser by 
making it propagate for a long distance under the influ- 
ence of a weak magnetic jlj or electromagnetic field 22 ■ 
This method has the advantage to enhance the phase 
shift by a long optical path without introducing costly 
strong target electromagnetic fields. In the case of a con- 
stant magnetic field on the order of 1 T, one encounters 
the limits of physical sensitivity to the QED nonlinear 
effects. In the case of an electromagnetic field, we may 
be sensitive to the QED-induced birefringence within a 
few days with a IJ CW laser according to the claim in 
[22'. However, if one aims at the sensitivity even beyond 
QED-induced birefringence as we discuss in section [3J 
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it is essential to introduce a high-intensity pulse even 
beyond the capability of the ELI facility [TW. In such 
circumstances the storage of a high-intensity laser pulse 
in a cavity is limited by the damage threshold of the 
optical elements needed to store the target field over a 
long time. 

On the other hand, if we could localize the field- 
induced refractive index change by tightly focusing a 
high-intensity target laser pulse and measuring the spa- 
tially inhomogeneous phase effect of the vacuum on a 
pulse-by-pulse basis, there will be no physical limit in 
increasing the intensity of the laser pulse until the vac- 
uum itself breaks down. In order to increase the shift of 
the refractive index, corresponding to the inverse of the 
phase velocities in Eq. ©, i.e., the intensity of the tar- 
get laser pulse as expected from Eq. ^ and Eq. (H)), it 
is necessary to use a focused laser pulse by confining the 
large laser energy into a small space-time volume. This 
causes a locally varying refractive index along the tra- 
jectory of the target laser pulse in vacuum. A variation 
of the refractive index arises over the high-intensity part 
and the remaining vacuum. If the probe laser penetrates 
into both parts simultaneously, the corresponding phase 
contrast should be embedded in the transverse profile 
of the same probe laser. Our suggestion is to directly 
measure the phase contrast and to determine the abso- 
lute value of the refractive index change by controlling 
the combination of polarizations of the probe and tar- 
get laser pulses. This should result in the birefringence 
as expected in Eq. The birefringence measurement 
based on the measurement of absolute phase velocities 
we suggest here should be contrasted to any existing 
techniques to measure the ellipticity where only relative 
phase differences can be discussed. 

In the following sections we introduce the basic con- 
cept of the phase-contrast Fourier imaging by crossing 
target and probe lasers and discuss a way to extract 
a physically induced phase in the presence of a phase 
background. We then discuss also other physical con- 
tributions beyond QED, to which this imaging method 
may be applied. 



2 Phase-contrast Fourier imaging 

We now consider an experimental setup, where we cre- 
ate a high-intensity spot by focusing a laser pulse in vac- 
uum and probe its refractive index shift by a second laser 
pulse. We call the first laser pulse the target laser pulse, 
while from hereon the second laser pulse will be denoted 
as the probe laser pulse. We need to detect the extremely 
small shift of the phase velocity by the target-probe in- 
teraction. For this we also need an intense probe laser 
in order to enhance the visibility. However, if we utilize 
conventional interferometer techniques, providing a ho- 
mogeneous phase contrast over the probe laser profile, 
such small refractive index changes are hard to detect. 
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Fig. 2 Far-field diffraction patterns from a thin-wire tar- 
get, when a Gaussian laser pulse is shot onto the wire. The 
left figure describes the situation without wire. The right fig- 
ure shows the case when the thin-wire target is horizontally 
arranged. 



This is because the resulting intensity modulation always 
appears on top of a huge pedestal intensity, with an ex- 
tremely small contrast between the modulation and the 
pedestal. Any photo-detection device will not be sensi- 
tive to the small number of photons spatially distributed 
over the pedestal intensity beyond 1 J (~ 10^^ visi- 
ble photons), due to the limited dynamic range of the 
photon intensity measurable by a camera pixel without 
causing saturation of the intensity measurement. On the 
other hand, broadening the dynamic range by lowering 
the gain of the electric amplification of photo-electrons 
degrades the sensitivity to the small number of the spa- 
tially distributed photons or the sensitivity to the small 
phase shift. Therefore, we need to invent a method that 
can spatially separate the weakly modulated character- 
istic intensity pattern from the strong pedestal. 

In order to overcome this difficulty, we suggest utiliz- 
ing the inhomogeneous phase-contrast Fourier imaging 
in the focal plane by focusing the probe laser. The phys- 
ically embedded phase contrast on the transverse pro- 
file of the probe laser amplitude is Fourier transformed 
onto the focal plane due to the effect of the added phase 
by the lens. Actually, a parabolic mirror is necessary to 
avoid dispersion and damage by high-intensity irradia- 
tion. This will be considered later. The intensity pat- 
tern in the focal plane exhibits the preferable feature, 
that the characteristic phase boundary causes outer re- 
gions of the intensity profile far from the focal point to 
expand, whereas a Gaussian laser beam with a homo- 
geneous phase converges into a small focal spot at its 
waist. It is instructive to illustrate the characteristic na- 
ture of the diffraction pattern from a wire-like target 
as shown in Fig. [21 Here the far-field pattern, known 
as Fraunhofer diffraction, is shown in the case when a 
Gaussian laser beam irradiates a thin- wire target. This 
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can be understood as the Fourier transform of the wire 
shape, approximated as a rectangle of 2fi x 2iy. It is well 
known that a lens produces a far-field diffraction pattern, 
corresponding to the exact Fourier-transformed image of 
the object in the front focal plane (e.g. see p51[?7] V In 
order to understand the diffraction image, we may qual- 
itatively refer to Babinet's principle, which states that 
the diffraction pattern from an opaque wire plus that 
of a slit of the same size and shape form an amplitude 
distribution identical to that of the incident wave. There- 
fore, the characteristic diffraction patterns from the wire 
and the slit are similar, but deviate from each other such 
that they interfere to reconstruct the incident wave. The 
intensity pattern after Fourier transform of such a rect- 
angular slit is expressed as 



definitions are as follows: 



(5) 



where lo^ 



^'f-y are the spatial frequen- 



Ts^ and - 
cies for the given position (x, y) in the focal plane of the 
lens/mirror with the focal length / at the wavelength 
A, respectively. In the case of a slit with fi ^ i', the 
rectangular profile in the focal plane emerges as a pat- 
tern of dark and bright fringes perpendicular to the slit 
(see Fig. [2] (right)). The narrower the slit size is, the 
further the fringes move apart. On the other hand, a 
Gaussian beam without wire or slit remains unchanged, 
because the Fourier transform of a Gaussian beam re- 
mains a Gaussian beam (see Fig. [5] (left)). This is the 
key feature that drastically improves the detectability of 
small phase shifts by sampling only outer parts of the 
diffraction pattern. This may also be interpreted as the 
counter-concept to the conventional spatial filter, where 
outer parts are eliminated to maintain a smooth phase 
on the transverse profile of the Gaussian distribution. 

Given the intuitive picture above, a quantitative for- 
mulation of our proposed method is presented as follows. 
In order to discuss the amount of the phase shift, we 
need a distinct geometry of both the target and probe 
lasers. Let us first consider the laser profile assuming 
Gaussian beams. The solution of the electromagnetic 
field propagation along z in vacuum is well-known [27 . 
The electric field component corresponding to the trans- 
verse mode l,m and e.g. polarized along y is expressed 
as E{x,y,z,t) = Re{eytpi,m{x,y, z)e^'^''} with 



'lpl^rn[X,y,Z) = A0——H1 I —^X I H,n \ —JT^y | X 



w{z) \ w{z) 
exp < —i[kz — {I + m + 1)77(0)] ~ ^'^ 



w{z) ' 

1 ik 
^ ^ 2R{z) 



where the Hi are l-th order Hermite polynomials, k = 
27r/A, r = ^/x"^ + J/^, wo is the waist, which cannot be 
smaller than A due to the diffraction limit, and other 



w{zf = wo'^ {1 + ] , 



zll 



ZR 



ri{z) = tan ^ ^ 



A 



(7) 
(8) 
(9) 
(10) 



In order to determine the normalization factor Aq, we 
use the orthonormal condition of the n*'' Hermite poly- 
nomial as follows 



+00 



With the replacement ^ = — a; we obtain 



(11) 



+00 



'V2 



-X e "0 dx / H 

Wo / J-00 



+ 00 



'V2 

Wo 



y]e ™o 



-u;o^7r2'+"Z!TO! (12) 

At z = we can then relate the amplitude Ao with the 
beam power Pq by 



Po 



yielding 



+00 



i'i.mix,y)'ipl^{x,y)dxdy (13) 



42 ^ ^ 

" ~ 27rcr22'+" 



l\m\ ' 



(14) 



where wq is replaced by 2a and Aq is the on-axis intensity 
at the waist. 

Figure [3] illustrates the conceptual experimental 
setup for the phase-contrast Fourier imaging. In what 
follows, subscripts p and t always denote the probe and 
target quantities, respectively. Both the target and probe 
laser beams are focused with different waist sizes Wto and 
with the incident beam diameters dt and dp, respec- 
tively. Both laser beams cross each other at their waists, 
where wave fronts are close to fiat with i? = 00 in Eq. ([5]). 
We assume that the target waist wot is smaller than the 
probe waist wop, which embeds the phase contrast S at 
zq within the amplitude on the transverse profile of the 
probe laser. The probe laser then propagates to a lens of 
focal length fp and the inverse Fourier imaging is per- 
formed in the back focal plane of that lens. 

We then define the geometry of the laser intersec- 
tion. Figure |4] illustrates geometrical relations where the 
tightly focused target pulse with time duration rj , beam 
waist diameter 2wot, and Rayleigh length Zfit propagates 
along the Z-axis, and the probe pulse with the larger 
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a) A< 2Wo, c) 




Fig. 3 Conceptual experimental setup for the suggested 
phase-contrast Fourier imaging. At the crossing point be- 
tween the probe and target laser beams, the target pulse 
causes a shift in the index of refraction, and a corresponding 
phase shift 5 embedded into the probe pulse as explained in 
Fig.H 

beam waist diameter 2w()p and longer time duration Tp 
propagates along the z-axis tilted by 9 with respect to 
the Z-axis. In FigUa) and b) the pulses are assumed to 
have a rectangular intensity profile along the propaga- 
tion direction. The dashed rectangular pulses occupy the 
positions at time t = and the solid ones are those at 
t — to. The probe and target pulses overlap each other 
at the position marked by *. At this moment in time, t 
is taken as zero. We need to express SI to estimate the 
pass length where an additional phase is embedded. The 
path length SI is defined as the distance where the front 
of the probe pulse meets the edges of the target laser 
at t ^ to, beyond which the target laser is no longer 
present. In FigU a) because A, B, and 61 form a right 
triangle, we hence obtain the following relation 



(15) 



where SI — cto, A ~ c(rf — to) and B — tan(7r— 0)c(rf — io 
with the velocity of light c, resulting in 

Tt 



to 



1 + cos(7r 



(16) 



Depending on the relation between the target beam 
waist 2'Wot and the pulse length ctj, namely, whether 
a) A < 2wot or b) where the probe wavefront meets the 
side of the target laser before reaching the tail as shown 
in FigHl the optical path length SI is expressed as 

a) 



SI — cto = 



cTt ^ 1 — cos 9 

for cTt < 2 Wot ^-5—7 (17) 



1 — cos 



b) 



SI ~ — — - tor CTt > iwot — 



COS0 



sm( 



smf 



(18) 



Fig. 4 Geometry of the embedded phase contrast in the 
probe pulse. 



where the equations should not be applied to the cases 
= or = TT. In the case a) the path lengths with 
phase shift are not constant below or above the star 
point along the wavefront of the probe laser. In the case 
b) the path lengths are constant over the part of the 
probe wavefront which penetrates both sides of the tar- 
get laser. The residual part along the wavefront, how- 
ever, meets the head or tail of the target laser pulse and 
causes deviations from the constant path length. An ex- 
actly equal path length over the probe wavefront during 
the propagation time of the target laser ~ 2z/i(/c is re- 
alized only in the case of 6* = 7r/2. In this case, after the 
penetration of the probe laser pulse, the profile of the 
probe laser in the x — y plane contains a trajectory with 
a constant phase shift S along the projection of the path 
of the target laser on the probe wavefront, as shown in 
Fig. He). 

We express the phase shift S in the vicinity of the 
waist \z\ — CTt 1 2 < znt, where we assume that the wave- 
front is fiat as indicated by Eq. ^ and Eq. ^ 



2-K 



SnSliptixp,yp), 



(19) 



where Sn is the refractive index shift, SI is the path 
length with an effectively constant phase shift over the 
crossing time and ipt{xp,yp) is a weighting function to 
reflect the path length difference depending on the inci- 
dent position with respect to the target profile expressed 
as a function of the position {xp,yp) in the transverse 
plane of the probe laser. If we limit the origin of the 
laser-induced refractive index change to QED, based on 
Eq. ([2]), Q, and dl]), we parametrize the refractive index 
shift as 



Sn 



qed 



C^o(l 



Et 



(20) 



where C is 4 or 7 for the polarization combinations || 
or _L, respectively. No is the coefficient to convert from 
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energy density to the refractive index shift defined as 
No = = 1.67 X 10-i2[;zm3/J]. The incident angle 

9 varies from to tt which is measured from the propa- 
gation direction of the target pulse to that of the probe 
pulse as depicted in Fig.H] Et is the energy of the target 
pulse given in [J], and niVQ^CTt is the volume in [/im^] for 
the given target profile with the waist wot from Eq. ([7]). 

By respecting the constant path length over ^ 2ziit/c 
for simplicity, we consider only the case of (jl8p with 
= tt/2. By substituting Eq. ^ and (UHl) into Eq. (fT9)) . 
we obtain the simplest expression for Sqed 



'^qed 



XpWotCTt 



(21) 



where CTt > 2wot must be satisfied from Eq. p8| and 
we take the approximation ipt{xp,yp) ^ 1 to simplify 
the following argument (if necessary, we may restore the 
target profile ipt {xp , yp) based on the precise profile of 
the target laser reflecting actual experimental setups). 
In this limit we approximate the target profile as a rect- 
angular of the size 2/i x 2v, inside which the phase shift 
is assigned to be constant. The effective slit sizes are de- 
fined by the transverse sizes of the focused laser beams 
through the relation 

jjL ^ Zf>t and V ^ wot- (22) 

We then explicitly define the window functions rec and 
rec as 



rec{ii, v) 
rec{^, v) 



1 for |a;| < and \y\ < v 

for |x| > [I or \y\ > v 

for \x\ < jj, and \y\ < v 

1 for \x\ > /i OY\y\> V 



(23) 



This window provides a unit region of a constant phase, 
which may be applied even to arbitrary phase maps com- 
posed of a collection of the unit window cells. 

We now discuss how the probe laser including the 
phase & embedded at the focal plane propagates into 
the image plane via the lens system and evaluate the 
expected intensity distribution in the image plane. Since 
our discussion is based on local phases with rectangular 
shape and both Fourier and inverse Fourier transforms 
of a rectangular function give identical sine functions, we 
represent the lens effect as Fourier transform. For each 
propagation from the object plane (xo,yo) at -^o to the 
image plane (x,?/) at z (see Figl3]), we always take the 
Fresnel diffraction. Based on © and the probe field 
profile in the plane where the phase b is embedded can 
be defined as 



T'(a;o,yo) = ApUi 



( xa 



V\/2cr 



( 2/0 



V\/2o 



(24) 



where Ap — 



— ,„,?^,, — r is the on-axis waist ampli- 
tude of the probe laser. The linearly synthesized ampli- 
tude at 00 is then expressed as 

<Z'(a;o, yo, ^o) = a(^ZQ)rec{[i, v)T{xq, yo) + 



(3{zo)rec{n,i^)T{xQ,yQ), (25) 

where a{z) and /3{z) are propagation factors of the probe 
waves at the point z after probe-target crossing. The 
functions a containing the phase shift 6 caused by the 
local refractive index shift and /3 are defined as 



a{z) 



- J{kz+5) 



(26) 



The Fourier transform F of the synthesized amplitude 
If" in the image plane {x,y) at z after the lens ^28j is 
expressed as 

F{^ixo,yo)}^ 

a{zo)F{rec{fi, v)T{xq, yo)} + /3{zQ)F{rec{fi, iy)T{xQ, yo)} 
= {a{zo) - l3{zo)) r rdxodyoT(xo,2/o)e-'("^""+"«^°) + 

J —fjj — 1/ 
nOO nOO 

(3{zo) / / dxodyoT{xo,yo)e-'^^-^''+'^-y'\(27) 



— OOJ —oo 



where we define (ujj;,ujy) = (-t^x, -M-y) at z. We in- 

jp^p jp^p 

troduce the coefhcient Cg for the first term in Eq. (P7)) . 
containing the information on how much the phase shift, 
representing the signal^ is localized, resulting in the 
photon-photon interaction. We decompose Cs into its 
real and imaginary parts, because Hermite polynomials 
contain even and odd functions and the non-zero values 
of these integrals appear even in the imaginary part. We 
denote them as 



ReCs(a;a;, Wy) = CsR 



dxQ / dyo{cos{uj.j:XQ) cos{ujyyQ) 



sin{u!xXo) sm{ujyya)}T{xo, yo) 



and 



dxo / dyo{cos{uJxXo) sm{u!yyo) + 
su\{uj.xXo) cos{uyyo)}T{xo, yo)- (28) 



We also define the coefficient Cb for the second term of 
Eq. (P7)) . which corresponds to the background pedestal 
as 

Cb{^x,^y) = 

{-if+-^{2a)^T:Hi{V2GUJx)Hm{V2<TUjy)e'i^ ,{2(i) 

where the fact is used that the nth-ovAev Hermite func- 
tion is the eigen-function of the Fourier transform, 
namely, {2Ti)-^/^F{H^{0e-^" /^} = {~iy^H,,{Lu)e'^" '\ 
We note that the coefficient (27r)~^/^ arises due to our 
definition of the Fourier transform with the prefactor of 
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unity applied to the lens system. The Fourier transform 
of the amplitude is then expressed as 

Fm - (a - I3){C,R + iCsi) + l3iCbR + iCti). (30) 

When no confusions are expected, we will omit 
{xo,yo,ZQ) and (u!x,ujy,z) for all relations below ((30)) . 
By substituting Eq. ^ and ^ into Eq. ((30)) . 

the intensity pattern at the image plane is expressed as 

\ Jp^P , 

{2(1 - cos ,5) (ac; - (CsrC^,b, + CiCm)) 
-2sin(5(C,/Cf,fl, - CsrCm) + C.C,*}, (31) 

where (/pAp)^^ arises from the Fresnel diffraction. Ac- 
cording to (HHj), when / + m is even or odd, Cm or C^r 
becomes zero, respectively. 

Equation ()31|) indicates that this method works as 
an interferometer via the cross terms with coefficients 
1 — cos (5 and sin^. The second term with sin (5 vanishes 
for any combinations of I and m, as long as the symmet- 
ric rectangular ranges around (a;o,J/o) = are assumed 
in the definitions of Eq. (j^Hj) . This interferometer differs 
from a conventional one in that the modulating part due 
to the phase shift b can be spatially separated from the 
confined strong part CbC^, due to the characteristic pat- 
tern of CgC* which creates images in the region of higher 
spatial frequencies. In Ea. ((3T)) the third term CtC^ cor- 
responds to the intense pedestal pattern insensitive to 
the phase S which keeps the same shape as that at zq 
with a different transverse scale based on Eg. ((^5)) . be- 
cause Hermite functions are the eigen-functions of the 
Fourier transform. The second term shows proportion- 
ality to S for (5 <C 1; however, the intensity pattern is 
constrained by Cb- This implies the phase information 
is attainable only in the vicinity of the background pat- 
tern CbC^, though the signal is strong. Therefore, the 
signal-to-pedestal ratio is not expected to be large. On 
the other hand, the first term indicates proportionality 
to (5^ for (5 <C 1 which implies a very weak signal; how- 
ever, the CgC* term is not affected by Cb and it produces 
a pattern characterized by spatial frequencies. Hence the 
signal-to-pedestal ratio is expected to be drastically im- 
proved circumventing the most intense spot. Therefore, 
depending on the value of 6 and the allowed dynamic 
range of the photo-detection device used, we have choices 
on which term we focus. 

We consider a Gaussian beam with I = m = 0. There- 
fore, we do not expect the terms proportional to phase 
S due to the imaginary parts originating from Hermite 
polynomials with odd orders. Nevertheless, if we need 
to stick to the proportionality to 6, we may add a local 
offset phase 7r/2 along the path of the target laser on 
the focal plane. We define this plane as the object plane 
for the inverse Fourier imaging as illustrated in FigjSj 
Such a setup can recover the sensitivity to the sign of 
the phase shift as well as the absolute value, because 
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Fig. 5 Conceptual experimental setup for the suggested 
phase-contrast Fourier imaging. At the crossing point be- 
tween the probe and target lasers, the target laser causes 
a shift in the index of refraction, which amounts to the re- 
fractive phase shift 5t,ac embedded into the probe laser as 
explained in Fig. |4| The offset phase is embedded by a holo- 
graphic phase plate in advance which provides the offset 
phase -|-7r/2 inside the rectangular region rec(/io, I'o) at zq 
via Fourier transform by the first lens. 



1 — cos((5 + 7r/2) in Eq. ()3ip becomes ^ 1 + 5. It is im- 
portant to be able to discuss whether the phase shift is 
increased or decreased, since it directly refiects the dy- 
namics of the local interaction. From a technical point 
of view, more importantly, this has a definite advantage 
of enhancing the signal due to the proportionality to b 
compared to the sensitivity in 1 — cos b in case of an 
extremely small b. However, in turn, one must accept the 
situation that the local offset phase contrast produces an 
intrinsic diffraction pattern as a new kind of pedestal, 
which now has an equal diffraction pattern compared 
to the one caused by the photon-photon interaction. 
Thanks to the proportionality to (5, we can reduce the 
intensity of the probe laser pulse. On the other hand, the 
new pedestal pattern would occupy the dynamic range 
of the camera device. In order to reduce the amount of 
the pedestal intensity, we may add more intelligent char- 
acteristic offset patterns by mixing 8 = b^ac i 7r/2 with 
the laser-induced vacuum phase shift i5„ac and the offset 
phase ±11/2 so that the offset diffraction patterns can 
destructively interfere in some points in the image plane 
thus still keeping a high signal-to-pedestal ratio. There- 
fore, the implementation of such local offset phases on 
the probe laser in advance is a key design issue, depend- 
ing on the dynamic range of the camera device. 

As illustrated in the zoom of the object plane (focal 
plane common to both lenses) in FigO we consider a 
rectangular phase offset so that it contains the region 
with phase b^ac ±7i'/2 in its center, that is, we define the 
offset region as rec(fiQ, vq) with fiQ — Nfi and vq = Nv 
and TV > 1 in Ea.((^ by giving the offset phase of ±7r/2. 
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Holographic device 




Fig. 6 An example how to implement the offset phase by a 
holographic device. 

For this region we introduce the coefficient Co(i^2:, Wj,) 
by replacing with (^Ojt'o) in Eg. ipS)) as well. The 

intensity profile in the focal plane is then re-expressed 
as 

{2{Co - Ct){Co ± S^acCs) + act}, (32) 

where ± are cases when the offset phase ±7r/2 are added 
in rec(/j,o, vq) and rec(/io, vq), respectively. We note that 
this relation is applicable to both real and imaginary co- 
efficients as long as either all real or all imaginary coeffi- 
cients are simultaneously zero. Actually we can confirm 
that Ea. ([5T|) becomes identical with Eg. ([5^ under this 
condition, when = 1 and Syac = 0, namely, Cg — Co 
and S = ±7r/2 are substituted into Ea.([3T|). 

The offset phase may be embedded by a sine distribu- 
tion which is, for example, produced via a step- like phase 
plate in advance via Fraunhofer diffraction by locating 
the plate at a far distance from the first lens in FiglS) 
This may provide the offset phase +tt/2 within the rect- 
angular region rec(/io, I'o) at zq via Fourier transform by 
the first lens in Figl5l If such a long distance is not avail- 
able, we may use a holographic device as illustrated in 
Figini where the phase plate is placed at the focal plane 
and a laser produces a proper Fourier image which is 
stored in the holographic device by mixing with a refer- 
ence laser. If we replace the reference laser by the probe 
laser in FigjSl we can produce the Fourier image in front 
of the first lens in FiglS] In a practical case as shown in 
FigM the holographic device may be located before the 
focal point when it records the local phases in advance, 
in order to supply an offset distance before the proper 
sine distribution is reconstructed by probe laser pulses 
at the exact point where we need it. 



3 Analysis in the Fourier image 

We performed numerical calculations with the rectangu- 
lar offset phase +7t/2 based on the setup illustrated in 
FigElwith Eq.dSll) for / m = (TEMqo). The param- 
eters used for Fig. [7] are summarized in Tab. [TJ where 
the parameters of the target and probe lasers, the em- 
bedded offset and the physical phase shifts due to the 
nonlinear QED effect used to obtain the Fourier trans- 
formed intensity distributions are specified. Figure [7] top- 
left illustrates the intensity pattern due to Eg. (15^ as a 
function (x, y) in the image plane when Syac — and 
the offset phase +n/2 is embedded only. The figure is 
plotted with an arbitrary unit for the contour height in 
logarithmic scale, by sampling values with 5 mm steps 
along the x and y-axes. Figure [7] top-right shows the 
expected number of pedestal photons, Np^d integrated 
over a 1cm x 1cm cell along the y-axis at a; = 0. Fig- 
ure [7] bottom-left shows Ngig — Nped per 1cm x 1cm cell 
along the y-axis at a; = 0, where N^ig is the integrated 
number of photons per 1cm x 1cm cell for signal, namely, 
with S = 6vac tt/2. In the actual experimental setup 
the subtraction should be performed on the shot-by-shot 
basis as illustrated in Fig|9] where a probe laser pulse is 
equally split into the signal path with the target laser 
pulse and the calibration path without it so that we can 
compare the two cases. The solid- red, dashed-blue and 
dotted-green histograms show the case when the probe 
wavelengths of 800nm, 840nm, and 760nm are assumed, 
respectively, and the yellow band shows the statistical 
fluctuations ^JNp^ due to the quantum efficiency of the 
photon detection. Figure [7] bottom-right shows the sta- 
tistical significance of the signal photons with respect 
to the statistical fiuctuations of the pedestal photons; 
{Nsig - Nped)/ ^jNped per 1cm x 1cm ceU along the y- 
axis at X = 0. The colors have the same meaning as those 
in Fig. [7] bottom-left. 

Figure [7] bottom-right indicates that we can expect 
several cells in which the number of signal photons is ei- 
ther increased or decreased by more than two standard 
deviations from the pedestal fiuctuations. If we could 
count the number of photons per 1cm x 1cm cell in the 
side band around the pedestal peak without detector sat- 
uration, we can measure the phase velocity shift due to 
the QED effect even by one probe-target laser crossing. 
This side band structure appears owing to interference 
between CQ{LOx,uJy) and Cs(wj;,Wy) in Eg.p2p. 

Although the spectral width of the probe laser shows 
faint effects as shown in Fig|7] bottom-left and bottom- 
right, the characteristic pattern along the y-axis is simi- 
lar. As long as the wavelength distribution can be mea- 
sured at the same time, we can reconstruct S^ac based on 
the measured wavelength distribution and the intensity 
pattern along the y-axis. 

The most difficult issue is the dynamic range of exist- 
ing cameras used in research which typically have 16-bit 
resolution and at most 28-bit per pixel. In order to solve 
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Fig. 7 Simulated distributions of V'o.oV'o.o (TEMoo) in image plane z based on Eq.((32]) with parameters given in Tab. [T] 
top-left: Patterns of Eq.((32]) with the offset phase +tt/2 only, top-right: The number of pedestal photons, Np^d integrated 
over a 1cm x 1cm cell along the y-axis at a; = of the top-left distribution, bottom-left: Nsig — Np^d per 1cm x 1cm cell 
along the y-axis at 1 = 0, where Nsig is the integrated number of photons per 1cm x 1cm cell with S = S^ac + Tr/2. The 
solid-red, dashed-blue, and dotted-green histograms in Fig[7] show the case when probe wavelengths of SOOnm, 840nm, and 
760nm are assumed, respectively, and the yellow band shows statistical fluctuations, -y/ Np^d, due to the quantum efficiency 
of the photon detector, bottom-right: The statistical significance of signal photons with respect to the statistical fluctuations 
of the background photons; {Nsig — Nped)/ \J Nped per 1cm x 1cm cell along the y-axis at a; = 0. The colors have the same 
meaning as those in the bottom-left chart. 
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Probe laser parameters 


n = 15 fs 


Tp = 22:Ht/c = 24 fs 


Et = 250J 


Ep = 25J 


At = 800 ± 40 nm 


Ap = 800 ± 40 nm 


dt — 39.8 cm 


dp — 7.0 cm 


ft = 75 cm 


fp = 25 cm 


wot ~ ^ = 0.96/im 




ZRt = TTTOoj/Af = 3.6/im 


^ijp = -KWop/Xp = 18.5^m 


Embedded physical phase by assuming only QED effect 


Sqed = 1.28 X 10"'" from Eq.lElJ with C = 4 and 6) = 7r/2 in Ea.(l20ll 


Shape of physical and offset phases 


fi = ZRt and f = wot 


Ho ~ 5 X jj, and uo = 5 x u with offset phase +n/2 



Table 1 Laser parameters used to produce Fig. [7] based on the conceptual experimental setup shown Figs. 2] and [S] The 
subscripts t and p refer to the target and probe lasers, respectively. This choice of parameters is explained in the text in 
sections [2] and [31 



the limited dynamic range, let us suppose that we sam- 
ple photons per 1cm x 1cm cell by ^ 10^ pixels. In such 
a case the number of photons per pixel is ^ 10^ with re- 
spect to ~ 10^^ photons at around 5cm from the pedestal 
peak (see Fig (7] top-right). Even if we use 10-bit reso- 
lution, the number of photons per resolution becomes 
10^/2^° - 10*^ photons. Compared to the N^tg - Nped 
of ~ 10^ at around 5cm from the pedestal peak (see 
FiglT] bottom-left) , the sensitivity of 10^ photons per res- 
olution is sufficient to observe the intensity modulations 
beyond two standard deviations from the pedestal fluc- 
tuations without intensity saturation (see FiglT] bottom- 
right). This suggests that in principle it is possible to 
detect the laser-induced QED effect by a single shot only 
if the conditions listed in TabUl are realized. Therefore, 
by assigning camera devices for individual 1cm x 1cm 
cell with ^10^ pixel readout, we can overcome the lim- 
ited dynamic range of cameras even with the currently 
existing technology. 

In order to study the laser-induced vacuum birefrin- 
gence, we inject a linearly polarized probe pulse whose 
electric field vector is turned by 45 deg with respect to 
that of the target pulse so that its electric field along the 
X and y axes are equal. We then put two polarization fil- 
ters at the image plane symmetrically with respect to 
y = as illustrated in FiglHlto cover the regions -\-y and 
—y along the y-axis, respectively, which select orthog- 
onal polarizations at the image plane. The asymmetry 
between the number of modulated photons from that of 
the pedestal pattern between the regions ±y provides 
direct information of the birefringence on the pulse-by- 
pulse basis. 

We note that this method bears similarity to that 
in [53], where two intense target laser pulses are treated 
as a matterless double slit and the interference between 
spherical waves from these slits is discussed as a signa- 
ture of the photon-photon interaction. In 23! the oc- 
currence of diffraction is caused by the laser-laser inter- 



action itself. In our method the target laser causes the 
refractive phase shift experienced by the probe laser, as 
indicated in Fig. l5] This phase shift is embedded in a 
refracted, nearly-plane wave in the forward direction of 
the probe laser, as explicitly formulated in Eq. (1^51) and 
Eq. (pS]) . We then set a lens to the right of the inter- 
action between the target and probe lasers as shown 
in Fig. l5l The diffraction or Fourier transform in our 
method is incurred by the added phase of the lens and 
the spherical wave propagation from the lens to the fo- 
cal plane. The advantage of our method is an enhanced 
sensitivity to a small phase shift on the pulse-by-pulse 
basis, as it is demonstrated due to a more efficient col- 
lection of photons by the lens using the much simpler 
target geometry. On the other hand, the disadvantage 
is the deviation from the ideal phases included in the 
path of the probe laser except the laser-induced vacuum 
phase. The ways to correct for this kind of background 
phase aberrations and the other background source for 
the phase-contrast Fourier imaging will be discussed in 
the following subsections. 

3. 1 Template analysis for local phase reconstruction 

In actual experiments it is unavoidable that the probe 
pulse includes local phase fluctuations on a pulse-by- 
pulse basis even in the absence of the laser-induced sig- 
nal 5 as illustrated in FiglHl The figure corresponds to 
the case when the phase contrast 5 in FiglSlc) is embed- 
ded in the presence of background phase fluctuations 
(pi = 4>{X) as a function of the position X = (xq, j/q) at 
the object plane where i denotes a corresponding region 
with the constant phase (pi in the transverse plane of the 
probe pulse. Compared to 5, the piS are expected to be 
much larger. However, if the values of the local phase set 
{<pi} on each probe pulse is suppressed below the offset 
phase 7r/2 and the set is a priori measured, we are in 
principle able to correct for the effect of the background 
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Fig. 8 The phase S induced by the target laser pulse in the 
presence of local phase fluctuations in the transverse plane 
of the probe laser pulse at the object plane where the target 
and the probe pulses cross each other. 

fluctuations. In the next subsection 13.21 we discuss how 
to measure the phase set on a pulse-by-pulse basis in 
detail. In this subsection, however, we focus on how to 
determine 5 on a pulse-by-pulse basis, if the measured 
phase set is given in advance. 

Let us extend the expressions from Eq. (l23l) through 
(P7|) . In general, the integration limits defined by 
Eq. (j23p and used in the first equation of (P7)) can take 
any shape and size. We replace the rectangular region rec 
with the region Ri = R{X), where a constant phase is 
mapped within Ri . By denoting the spatial frequency as 
W = {uj^,ujy) = {2nx/{fp\p),2Try/{fpXp)), for the po- 
sition (x, y) at the image plane with the integral kernel 
f{W,X) = r(a;o,2/o)e"'^"^'^°+'^"^''\ the Fourier trans- 
form including the local phase fluctuations is ex- 
pressed as 



the laser-induced phase is given by 

hgiW; (/)) = {Re^iW; 0)}^ + {Im^{W; 0)}^. (34) 

We now include as well a template of the laser- 
induced phase by the target laser pulse 6 = 6{X). The 
phase shift S can be evaluated from the geometry of the 
energy density profile of the target laser pulse. Based on 
Eq. p^ we parametrize S as 

5 = taptiX), (35) 

where k is a constant parameter that considers the abso- 
lute value of the phase shift induced by the target laser. 
The profile can be a priori determined by the experi- 
mental design of the focal spot. We can monitor if the 
center of the spot is in fact stable and further correct for 
its deviation from the fixed geometry of the target laser. 
Given 6, we only have to replace the phase by 0i — (pi+S 
with a constant parameter k as follows 

Ibg+sig{W;(j) + Kifit) = 

{Re<f (T^; + Kipt)V + {lm<Ir{W; 4> + Kipt)f, (36) 

where bg + sig refers to the fact that the laser-induced 
phase is embedded in the background phase fiuctuations. 

Given the measured intensity pattern Imeas in the 
image plane per probe pulse, we define with Eq. 
as a function of k 

1 \Irneas{Wj) - hg+sigjWj-, (j> + KLpt)\'^ 

Nw-l^ I-meas{Wj)+hg+sig{Wf,(l) + Kipt) 

where Nw is the number of sampling points in the im- 
age plane and j runs over all regions in this plane. The 
parameter k can be determined by minimizing on a 
pulse-by-pulse basis within the required accuracy. 



^{W;ct))^F{^{X-cty)) 

Nx „ 

= 5^{a(0,) - /?} / dXf{X,W) 

Nx 

= 5^{a(<^0 - I3}T,{W) + l3Too{W), (33) 

i 

where Nx is the number of regions in the transverse 
plane at z, a{<j),) = e'(fc2+</'.)^ p ^ e'*=^ Ti{W) = 
J^^dXf{X,W), and TooiW) = J^^dX f{X,W). We 
note that this expression corresponds to the regional cut 
and paste on T{xq, yo); i.e., cutting a region with a phase 
determined from /? at z and paste the same region by 
adding 4>i in a{(j)i). 

Given on a pulse-by-pulse basis, we can numeri- 
cally calculate the real and imaginary parts of !f'(PF; 0). 
The estimated background intensity pattern Ibg{(j}) in 
the image plane with the phase fluctuations 4> without 



3.2 Corrections on pulse-by-pulse phase aberrations 

Figure OiUustrates a schematic view of the entire system 
for the phase-contrast Fourier imaging including parts to 
correct all phase aberrations in the system. The target 
laser pulse moves perpendicular to the drawing plane. 
Its focus or waist lies in this plane. The signal path 
(SP) consists of the inverse Fourier transform part as 
discussed in FigIS] and an array of mega-pixel camera 
sensors at the end to sample the intensity profile by 
individual 1cm x 1cm cells as discussed with Fig|7l In 
the SP the probe laser pulses are injected with the po- 
larization tilted by 45 deg with respect to that of the 
target laser pulses. In front of the sensors, two polar- 
izers (P) selecting photons with the orthogonal combi- 
nation of polarizations so that the birefringence can be 
measured in a shot, which allows the statistical integra- 
tion of the measurement over many shots by minimiz- 
ing the systematic error due to shot-by-shot fluctuations 
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of the probe pulse energy. After the implementation of 
the holographic plate (HP) which produces the offset 
phase contrast in the laser interaction zone, we intro- 
duce a beam splitter (BS2) followed by the identical im- 
age transferring system as that in the SP. We refer to this 
leg as the calibration path (CP). We classify the origins 
of local phase fluctuations into the static component by 
the optical elements in the paths and the pulse-by-pulse 
component such as wavefront fluctuations included in 
the probe pulse coming from the upstream laser system. 
Since the repetition rate of the target laser is limited, 
we may inject a single-mode CW laser with the same 
wavelength as the dominant part of the probe pulse spec- 
trum into both the SP and the CP, while the target laser 
pulses are not injected (Wavefront aberrations resulting 
from reflection by BSl are predetermined). Because a 
weak CW laser may be realized as a perfect Gaussian 
beam running in the TEMqo mode at a single longitu- 
dinal mode, we expect to be able to accurately deter- 
mine the static phase component as the average value 
by using huge photon statistics accumulated over a long 
time period for an experiment while probe pulses are not 
injected. For the pulse-by-pulse component we use the 
intensity profile observed at the end of the CP to recon- 
struct a set of local phases caused by wavefront fluctu- 
ations included in the probe pulse on the pulse-by-pulse 
basis. 

The measurable four types of local phase sets are 
denoted as (j)g^, 'f^cp^ ^'^'^ '^cp'^ where super- 

scripts specify cases of CW and pulse laser injections, 
respectively, and subscripts refer to the different paths 
the beams take. In the following, all phase sets are in- 
terpreted as those defined on the focal plane, even if the 
local phases are actually embedded in different prop- 
agation points. The two phase sets in the SP are ex- 
pressed by phases (/s's with subscripts corresponding to 
the names of the optical elements along the path in FiglH] 
as follows: 

05]f = ipHP + VBS2 + VPMIsp + VPM2sp + VPspX^^) 

and 



,PLS 
i^SP 



IPPLS 



,CW 
l^SP > 



(39) 



where BSl should be removed when the SP is active, 
and (fipLS is the pure phase set caused by only the pulse- 
by-pulse component which is not correctable by the CW 
laser. The two phase sets in the CP are expressed as well 

•^CP = VHP + VBS2 + VPMlcP + VPM2CP + V'Pcp '(40) 

and 



This implies that (t)gp^ can be restored by other measur- 
able quantities, which is a necessary condition to allow 
the correction within the same probe pulse injection in 
the SP in the presence of the laser-induced vacuum phase 
shift. We finally describe the entire phase set in the focal 
plane in the SP when a target laser pulse exists as 



(43) 



where Syac = Kft as parametrized in Ea. (|35p . By sub- 
stituting Ea. (|43p into Eq. p7p . we can, in principle, de- 
termine K for the physical template based on the target 
laser profile tpt- 

The template analysis discussed in Sect. I3.1l can also 
be applied to determine the individual set of phases in 
the right hand side of Ea. (|l^ . By assigning a square 
shape to the region Ri in Eq. (j33p , representing a cell in- 
stead of physical template (ft in Eq. ([ST)) , we estimate Ki 
for each Ri{X). The number of photons at a point Wi in 
the image plane contains the convoluted phase informa- 
tion of the amplitude from all points in the transverse 
plane of the probe Nx as seen from Eq. ([55]) . Therefore, 
as long as the number of sampling points in the image 
plane Nw is larger than that in the transverse probe 
profile iVx, we can, in principle, determine a phase set 
from Eq. (j37p by scanning Ki over the expected dynamic 
range of the phase variation. The achievable resolution 
of the phase reconstruction depends on the scanning step 
on Ki in the x^-test. As discussed with Figl7]the phase- 
contrast Fourier imaging achieves at least the sensitiv- 
ity of ~ 10""'^° for the physical phase shift by sampling 
the side band of the intensity distribution on the image 
plane. Therefore, we can introduce the same resolution 
step to determine Ki. We may measure the initial coarse 
phase sets a priori by a commercially available wavefront 
sensor. From the phase measurement we can extract the 
set of phases at the focal plane by performing Fourier 
transform from the image plane back to the focal plane. 
Starting from this initial phase set at the focal plane, 
we perform the x^-test to determine Ki more accurately 
by comparing the computed Fourier image to the mea- 
sured intensity at the image plane. If the resolution of 
commercially available wavefront sensors is limited to 
^ A/100, we would need to repeat the two-dimensional 
inverse Fourier transform from the focal plane to the im- 
age plane more than 10* times for scanning Ki, in order 
to reach the same phase resolution as ~ 10"^". Accord- 
ingly, a proper computing power is necessary to restore 
the sets of the offset phases in Eq. l|l^ on the pulse-by- 
pulse basis. 



iPLS 
l^CP 



Vpls ■ 



• CW 



(41) 

3.3 Background in the phase- contrast Fourier imaging 



Combining Eg. ((39)) and ((4T|) . we can restore the offset 
phase for the probe pulse injection in the SP, (j)sp^ by 
the other measured sets of phases as 



,PLS 
l^SP 



,PLS 
''CP 



,cw 

''CP 



• CW 
^SP ■ 



(42) 



A background source of the current measurement is the 
refractive index shift due to the plasma creation from the 
residual gas along the path of the focused target laser 
pulse. The refractive index of the static plasma in the 
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Fig. 9 Setup to correct local phase fluctuations. 

limit of negligible collisions between charged particles is 
expressed as 

iV = Jl-^, (44) 

where wq is the angular frequency of the target laser, ujp 
is the plasma angular frequency defined as ■\/47re^ne/me 
and 7 is the relativistic Lorentz factor given as y^l + Qq 
with ao = 0.85 x IQ-^ X[fim] a/Iq [W/cm'^] . In the low- 
pressure limit of the residual gas, the amount of refrac- 
tive index shift AN = iV — 1 is expressed as ujp'^ /2jujq^ . 
Although the refractive index in the plasma becomes 
smaller than that of the peripheral area with neutral 
atoms, the inverted phase contrast of the phase shift in- 
side the probe pulse still maintains a rectangular shape 
along the trajectory of the target laser. Therefore, it 
should produce the characteristic diffraction pattern at 
similar locations to the nonlinear QED case as expected 
from the Babinet's principle, which requires that the 
diffraction pattern from an opaque slit plus the inverted 
slit of the same size and shape form an amplitude dis- 
tribution identical to that of the incident wave as we 
discussed in section In order to reduce this effect, we 
need to reduce the electron density in the residual gas. 
If we take 7 '-^ 1 as the upper limit of the AN estimate, 
the refractive index shift ~ 10~^^ due to the nonlinear 
QED effect for a reference energy density ^ IJ/zim'^, cor- 
responding to a residual gas pressure of ~ 10~^ Pa. The 
collisional frequency due to interactions between elec- 
trons and ions is expected to be 10^— lO^s^^ at the crit- 
ical electron density ncr[cm~^] = 1.12 x 10^^/A^[/im], 
where LOp equals wq. For a duration time of ~fs of the 
target laser pulse, the inverse bremsstrahlung radiation 
due to collisional processes in the residual gas is negligi- 
ble at - 10"^ Pa. 

Plasma formation is also caused by the probe pulse 
along its waist over a distance of ^ 1 mm. At a pressure 
of ~ 10^^ Pa the associated plasma induced phase shift 



is one order of magnitude smaller than that due to QED. 
Hence the pressure of the residual gas in the interaction 
chamber has to be kept at this level. 

We note that the actual processes will be more dy- 
namical, due to the pondermotive force executed by the 
high-intensity laser field. In such a case the refractive 
index shift based on static plasma gives only the upper 
bound on the amount of the local refractive and phase 
shift. 

4 Potential efTects beyond QED 

In the previous sections we discussed the design of 
the phase-contrast Fourier imaging by aiming at prob- 
ing the vacuum birefringence through the QED effect, 
namely, the electron-positron loop to which photons cou- 
ple. However, if the quark mass in vacuum is of the same 
order as the electron mass, we should expect that quarks 
also contribute to the vacuum birefringence by replac- 
ing the electron-positron loop with the quark-antiquark 
loop. Whether this effect has a sizable contribution or 
not is, however, difficult to quantify with presently ex- 
isting field theoretical approaches, because of the strong 
coupling of quantum-chromodynamics (QCD) in vac- 
uum, where the coupling is too large to allow for a per- 
turbative treatment. Moreover, bare quark masses not 
confined in hadrons are not precisely known. In addi- 
tion to calculations based on the QCD field theory [16l 
[T7] , there is another possibility that the duality between 
string theory with higher dimensions and field theory 
in 3+1 dimensions (holography) [TBlfT^ could be di- 
rectly applicable to this birefringence problem [JO]- The 
QCD and holographic approaches may give different pre- 
dictions for the balance of coefficients between the two 
terms of the Euler-Heisenberg Lagrangian. Therefore, we 
may be able to pin down such theoretical issues by ac- 
cumulating statistics more than a single shot and also 
expecting further increase of the laser intensity in the 
future. 

Moreover, we note that because the photon-photon 
scattering cross section of QED interaction in the per- 
turbative regime is so small, 10^*'^ b at optical frequency 
(see [29ll30] ). we experience little 'noise', providing a 
pristine experimental environment to search for some- 
thing beyond QED. Suppose then the detected disper- 
sion and birefringence quantitatively deviate from the 
expectation of QED, including potential QCD correc- 
tions. This should indicate that undiscovered fields may 
be mediating photons beyond QED and QCD. Scalar 
and pseudoscalar types of fields in vacuum may con- 
tribute via the first and second products in the brackets 
of Eq. ([T]) , respectively. They may be candidates of cold 
dark matter, if the coupling to photons and the mass 
are reasonably small |31) . Therefore, the measurement 
of the absolute value of the phase shift depending on the 
polarization combinations and the comparison to the ex- 
pectations from nonlinear QED including potential QCD 
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corrections may be a general test of unknown nature in 
vacuum. 



5 Conclusion 

We suggest an approach to probe the vacuum birefrin- 
gence under the influence of intense lasers. The phase- 
contrast Fourier imaging technique can provide a sensi- 
tive method to measure the absolute phase shift of light 
crossing intense laser fields. With this method nonlinear 
QED effects of the Eulcr-Hciscnberg Lagrangian may be 
detected requiring no more than lasers of the hundred 
PW-class. The method provides a window for scoping 
the vacuum via the dynamics of the electron mass scale 
and possibly the lightest quark mass. Such a detection 
has never been made to date, and it heralds the research 
in the physics of the vacuum with a high-field approach. 
Given the high-intense optical lasers available in the ELI 
project |15| in the near future, the realization of this 
suggestion may become an exciting challenge for future 
experiments exploring vacuum physics. 
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